typeset using JPSJ.sty <ver.l.Ob> 



Random Fixed Point of Three-Dimensional Random-Bond Ising Models 



Koji Hukushima 



o 
o 
o 

(N 



Institute for Solid State Physics, Univ. of Tokyo, 7-22-1 Roppongi, Minato-ku, Tokyo 106-8666 

(Received November 12, 1999) 

The fixed-point structure of three-dimensional bond-disordered Ising models is investigated 
using the numerical domain- wall renormalization-group method. It is found that, in the ±J Ising 
model, there exists a non-trivial fixed point along the phase boundary between the paramagnetic 
and ferromagnetic phases. The fixed-point Hamiltonian of the ±J model numerically coincides 
with that of the unfrustrated random Ising models, strongly suggesting that both belong to the 
same universality class. Another fixed point corresponding to the multicritical point is also found 
in the ±J model. Critical properties associated with the fixed point are qualitatively consistent 
with theoretical predictions. 
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The influence of quenched disorder on model systems 
has attracted considerable interest in the field of statis- 
tical physics. The first remarkable criterion was given 
by Harris J$ who claimed that if the specific- heat expo- 
nent Qfpurc of the pure system is positive, disorder be- 
comes relevant, implying that a new random fixed point 
governs critical phenomena of the random system. One 
of the simplest models belonging to such a class is the 
three-dimensional {3D) Ising model. Critical exponents 
associated with the random fixed point have been iowas.- 
tigated for dilution- type disorder by experimentalpLj'tf 
theoreticafiti and numerical approaches .LrcP Recently, 
extensive Monte Carlo (MC) studied^ have clarified that 
the critical exponents of the 3D site-diluted Ising model 
are independent of the concentration of the site dilu- 
tion p, suggesting the existence of a random fixed point. 
These results were obtained by carefully taking into ac- 
count correction for finite-size [Scaling, unless the ex- 
ponents clearly depended on p.a 1 Experimentally, the 
critical exponents of randomly diluted antiferromagnetic 
Ising compounds, Fe :E Zni_ :E F2EiQ' and Mn^Ziii-aFaP 
arc distinct from those of the pure 3D Ising model. 

While the existence of the random fixed point has 
been established for the 3D site-diluted Ising model, 
the idea of the universality class for random systems, 
namely classification by fixed points, has not been ex- 
plored yet as compared with various pure systems. In 
particular, the question as to whether the random fixed 
point is universal irrespective of the type of disorder or 
not is a non-trivial problem. In the present work, we 
study critical phenomena associated with the ferromag- 
netic phase transition in 3D site- and bond-diluted and 
±J Ising models. The main purpose is to determine 
the fixed-point structure of 31? random-bond Ising mod- 
els by making use of a numerical renormalization-group 
(RG) analysis. Our strategy is based on the domain-wall 
RG (DWRG) method proposed by McMillan H-EJ This 
method has been ..applied to a 2D frustrated random- 
bond Ising modelBH where there is no random fixed 
point, and recently to a 2D ± J frustrated random-bond 



three-state Potts model0' which displays a non-trivial 
random fixed point. In this paper, we show systematic 
RG flow diagrams for 3D random Ising spin systems, 
which convinces us of the existence of the random fixed 
point. 

Let us first explain briefly the RG scheme rfpd the 
numerical method used by us. In the DWRG,BEP the 
following domain-wall free energy AFj of a spin system 
on a cube with the size L is regarded as an effective 
coupling associated with a length scale L for a particular 
bond configuration, denoted by J; 

« = h M (i) 

where Zp(ap)(F) is the partition function of the cube 
at temperature T under (anti-) periodic boundary con- 
ditions for a given direction, while the periodic bound- 
ary conditions are imposed for the remaining directions. 
In disordered systems, the distribution P(AF) of the 
effective couplings over the bond configurations is con- 
sidered to be a relevant quantity. Therefore, in the 
DWRG scheme, we are interested in how the distribu- 
tion is renormalized as L increases. In the ferromag- 
netic phase, the expectation value of the distribution ap- 
proaches infinity as L increases, while it vanishes in the 
paramagnetic phase. Fixed points are characterized by 
an invariant distribution of the coupling under a DWRG 
transformation, namely, increasing L. For example, the 
unstable fixed-point distribution corresponding to the 
pure-ferromagnetic phase transition is a delta function 
with a non-zero mean. The random fixed point, if any, is 
expected to have a non-trivial distribution with a finite 
width. 

It is convenient to consider typical quantities charac- 
terizing the distribution P(AF), instead of the distribu- 
tion itself. We discuss here the mean AF and the width 
<t(AF) of the distribution which are evaluated as 



AF = / dAF AF F(AF), 



(2) 
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a 2 (AF) = AF 2 — AF 



(3) 



respectively. Using these quantities, we define two re- 
duced parameters, ■^.~aLAF)/AF and t = T/ AF, as in 
the previous works B 113 liJ The renormalized parameters 
(r(L),t(L)) of the length scale L are estimated numeri- 
cally with bare parameters in a model Hamiltonian fixed. 
Then, the RG flow is represented by an arrow connected 
from the point (r(L),t(L)) of size L to (r(L'),t(L')) of a 
larger size L' . Fixed points should be observed as points 
where the position is invariant under the RG transfor- 
mation. Linearizing the flow about a fixed point (r* , £*), 
we obtain 



r{L')-r* 
t(L') - t* 



= T 



r(L) - r* 
t{L) - 1* 



(4) 



where T is a RG transformation matrix whose eigen- 
vectors are scaling axes of the RG flow. The matrix 
elements of T as well as the fixed point (r*,i*) can be 
determined by least-squares fitting from the data point 
numerically obtained. The critical exponents y are ob- 
tained as log A/ log(L'/L) with the eigenvalue A of T. 

We consider Ising models with quenched disorder 
which is defined on a simple cubic lattice. The model 
Hamiltonian is 



(io) 



(5) 



where the sum runs over nearest-neighbor sites. In order 
to obtain the domain-wall free energy, we use a recently 
proposed MC methocO which enables us to evaluate 
the free-energy difference directly with sufficient accu- 
racy. In the novel MC algorithm, we introduce a dynam- 
ical variable specifying the boundary conditions, which is 
the sign of the interactions between the first and the last 
layer of the cube for the given directipa. The algorithm is 
based on the exchange MC method Jl^ sometimes called 
the parallel tempering ,ES which turns out to be reason- 
ably efficient for randomly frustrated spin systems. 
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Fig. 1. Flow diagram for the 3D unfrustrated bond-diluted 
Ising model. The random fixed point is found numerically at 
(0.63(2), 1.77(5)). The bold arrows indicate the eigenvectors, 
whose exponents are y± = 1.47(4) and yi = —1.3(4). 

First we study a 3D bond-diluted Ising model, which 
is given by eq. (||) with the bond distribution, P(Jij) — 
pS(Jij — J) + (1 — p)S(Jij). The bond concentrations 



studied are p = 1.00, 0.90, 0.65, 0.45 and 0.35 with the 
system sizes L = 8 and L — 12. Sample averages are 
taken over 128 — 1984 samples depending on the size and 
the concentration. Errors are estimated from statistical 
fluctuation over samples. The number of temperature 
points in the exchange MC is fixed at 20. We distribute 
these temperatures to replicas for each p such that the 
acceptance ratio for each exchange process becomes con- 
stant. We show the RG flow diagram in Fig. 1. When 
disorder is absent, corresponding to the x — axis in 
Fig. 1, the pure fixed point is observed at (0, 1.63), de- 
noted by P. Near the fixed point P, the arrows flow 
away from P as disorder is introduced. This means that 
the pure fixed point of the 3D Ising model is unstable 
against the disorder, consistent with the Harris crite- 
rionlP Meanwhile, this finding suggests that a charac- 
teristic feature of the RG flow is reproduced within the 
system sizes studied. Apart from the pure fixed point 
P, we find another fixed point, denoted by R, along the 
phase boundary. The RG flows approach R from both 
sides, indicating that it is an attractive fixed point along 
the critical surface. The fixed point R governs critical 
phenomena of the disordered ferromagnetic phase tran- 
sition. Therefore, it should be called the random fixed 
point. The position of R in the parameter space is ob- 
tained numerically as (0.63(2), 1.77(5)). At R, the expo- 
nents yi = 1.47(4), y2 = —1.3(4) and the corresponding 
eigenvectors are indicated by the bold arrows in Fig. 1. 
The critical exponent the inverse of the larger eigen- 
value ux , is compatible with that of the site-diluted Ising 
modelja) although our estimate is not so accurate. All 
the arrows below T c go to the unique ferromagnetic fixed 
point at (0, 0). 
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Fig. 2. Phase diagram of the 3D random-bond Ising models with 
the bond (triangle) and the site (box) dilution. In the inset, a 
schematic flow diagram for the diluted Ising models is shown as 
a function of temperature and disorder. The fixed points and 
the RG flow are indicated. 

The transition temperatures_are estimated by a naive 
finite-size scaling assumption J§t°P AF = f((T-T c )L 1 / v ) 
with the correlation-length exponent namely, the 
crossing point of AD with L — 8 and 12 as a function of 
temperature is located on T c . This scaling form is similar 
to that of the Binder parameter frequently used. The es- 
timated T c is shown in Fig. 2. The exponent v obtained 
by the scaling depends significantly on the concentration 
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p, similar to those observed as an effective exponent in 
the MC simulation of the site-disordered Ising model© 
It is found from the RG flow diagram that the system 
has a subleading scaling parameter which gives rise to a 
systematic correction to the scaling. Therefore, we con- 
clude that the continuously varying exponent observed 
in the naive scaling is due to the subleading parameter. 
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Fig. 3. Flow diagram for the 3D unfrustrated site-diluted Ising 
model. The random fixed point is located on the point 
(0.64(4), 1.71(9)), whose exponents are y\ = 1.37(9) and y2 = 
-0.9(5). 

We also investigate a 3D site-diluted Ising model by 
the same procedure as in the bond-diluted Ising model 
described above. The model is given by eq. (||) but 
with the bond Jij = Jeitj and the distribution, P{ti) = 
pS(ei — 1) + (1 — p)6(ei). As shown in Fig. 3, there exists 
a random fixed point, which is consistent with the uni- 
versality scenario observed in this 3D site-diluted Ising 
model© In the MC simulation,!!^ the correction to the 
scaling becomes smaller around p — 0.80. This can be 
explained by the finding that the random fixed point 
we found is located on the position near the bare cou- 
pling with p = 0.80. An interesting point to note is that 
the position of the fixed point is close to that observed 
in the 3D bond-diluted Ising model, though the corre- 
sponding bare parameters such as the concentration p 
and the temperature differ between these two models, as 
seen in Fig. 2. This finding implies that the random fixed 
point is universal for a large class of the 3D unfrustrated 
random Ising models. 

Next we consider a 31? ± J Ising model, where the in- 
teractions Jij are randomly distributed according to the 
bimodal distribution, P(Jij) = p5(Jij — J) + (l—p)5(Jij + 
J). The multi-spin coding technique can be easily im- 
plemented in this model. For that purpose, we consider 
32 temperature points in the exchange MC simulation. 
In this model, the ferromagnetically ordered state sur- 
vives at low temperatures up to a critical concentration, 
recently estimated at p c — 0. 7673(3), B while below p c a 
spin glass (SG) phase appears. We estimate the domain- 
wall free energy with L — 8 and 12 for a wide range 
of concentration p including the critical concentration. 
Sample averages are taken over 240 — 6800 samples de- 
pending on the size and the concentration. 

We show in Fig. 4 the RG flow diagram for the 3D ± J 
Ising SG model. The arrow connects results for L = 8 
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Fig. 4. RG flow diagram for the 3D± J Ising SG model. The pure 
critical fixed point is denoted by P, the random fixed point by R 
and the multicritical fixed point by TV. The bold arrows repre- 
sent the eigenvector of the RG transformation matrix T at R and 
N . The broken line represents a critical surface as a guide for the 
eyes. The random fixed point R at (0.66(2), 1.74(1)) is charac- 
terized by the eigenvalues y\ = 1.52(2) and y2 = —0.42(13). The 
exponents associated with the fixed point N at (1.66(6), 1.17(2)) 
are y\ = 1.2(2) and y2 = 0.62(6). Data points on the Nishimori 
line are marked by gray symbols. 

and L' = 12. Here, we again find the random fixed 
point R as an attractor along the ferromagnetic phase 
boundary. The position of R is also close to that ob- 
served in the diluted Ising models. This finding indi- 
cates that both belong to the same universality class as 
the ferromagnetic phase transition. In other words, the 
renormalized coupling constants (or coarse-grained spin 
configurations) are independent of the details of the mi- 
croscopic Hamiltonian and of whether it is frustrated or 
not. One of the inherent characteristics of the ± J model 
is the existence of a highly symmetric line, which we 
call the Nishimori line.EEP There were ^PYPral theoretical 
works concerning the Nishimori lincMBBBEP It was 
suggested by the e-expansion method and a symmetry ar- 
gumentLjtj' that the multicritical point must be located 
on the line. This was confirmed by MC suniilationsla'cll' 
and high-temperature series expansions .USES As shown 
in Fig. 4, we also observe a fixed point, denoted by N, 
corresponding to the multicritical point, where both scal- 
ing axes have a positive eigenvalue. One of the scaling 
axes governed by the larger eigenvalue at N almost coin- 
cides with the Nishimori line, whichJ-s-irj agreement with 
the predictions by the e-expansion. t3[3 

As seen in Fig. 5, the estimated critical tempera- 
tures for the concentrations simulated are consistent with 
those obtained-,by. the large-scale MC simulations up to 
size L = 101.ESE3) Our analysis can also be performed 
for the SG transition at smaller values of p. The SG tran- 
sition temperature Tgc is determined from the crossing 
of <t(AF). The value of T SG at p =p0.50 is in good 
agreement with the recent estimation. c3> Along the fer- 
romagnetic phase boundary, it is natural to expect that 
there exists another fixed point Z at zero temperature, 
which separates the ferromagnetic phase from the SG 
one. Because the eigenvalues of the RG matrix at N are 
positive for both scaling axes, the thermal axis would be 
irrelevant in contrast to the diluted Ising models, namely, 
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an irrelevant flow into the zero-temperature fixed point 
Z originates at N, as shown in the inset of Fig. 5. The 
zero-temperature fixed point governs the critical behav- 
ior between the ferromagnetic and SG phase. In the 
present study, however, we cannot reach temperatures 
that are sufficiently low to extract the zero-temperature 
fixed point. A modern optimization technique such as 
a genetic algorithm would be useful in searching for the 
location of the fixed point. 
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Fig. 5. Phase diagram of the 3D ± J Ising SG model. Open tri- 
angles are transition temperature estimated in the present work. 
The SG transition temperatures marked by the solid triangles at 
p = 0.50 and 0.76 are determined from <j(AF). The open cir- 
cles reppme«t estimations by the nonequilibrium MC relaxation 
method|i3'E2P • The dotted line represents the so-called Nishimori 
line. A schematic flow diagram for the ± J model is drawn in the 
inset. The fixed points R and N, and the RG flow around them 
are supported in the present work. 

Now we comment on the correction to the finite-size 
scaling. In order to obtain the critical exponent, one 
uses a range of the system size frequently independent 
of the bare parameters in the model Hamiltonian. When 
the system has a subleading scaling parameter, it causes 
systematic corrections to the scaling. However, it is not 
known how such corrections affect the leading scaling a 
priori. The RG flow diagram gives a good indication of 
necessity and justification for corrections to the leading 
scaling term. Once the fixed-point structure of inter- 
est is clarified from the flow, one of the best numerical 
approaches to obtain the critical exponents, within re- 
stricted computer facilities, is to choose a model param- 
eter close to the fixed Hamiltonian. 

In experiments, the 3-D ± J Ising model is approx- 
imately realized in a mixed compound with strong 
anisotropy such as Fe^Mni-^TiOa. Our findings sug- 
gest that the critical phenomena of the compound 
Fe^Mni^TiOs near x ~ I belong to the same univer- 
sality class as the 3-D random Ising models. We expect 
that various random Ising systems such as a mixed Ising 
spin compound Fe a ,Coi_ 2 ,F2,ES' though beyond the scope 
of he present study, are also categorized into the same 
universality class. 

In conclusion, we have studied the fixed-point struc- 
ture of the 3D random-bond Ising models using the nu- 
merical DWRG method. We have found in the ± J Ising 
model, a random fixed point besides the pure and mul- 
ticritical ones along the ferromagnetic phase boundary. 



Furthermore, the observed random fixed point is found 
to be very close to that of the 3D non-frustrated dilute 
Ising models, while bare parameters in the model Hamil- 
tonians differ entirely. This fact strongly suggests that 
there exists a universal fixed point characterizing the 3D 
disordered ferromagnetic Ising model irrespective of the 
type of disorder. The present work is, to our knowledge, 
the first of its kind performed for determining the fixed- 
point structure of a 3D random spin system. We consider 
that the present numerical RG approach is quite useful 
for understanding the universality class of random spin 
systems, including spin glasses. 
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